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<1> . Abstract 

Deep-inelastic scattering data, taken with the HI detector at HERA, are used to study the 
event shape variables thrust, jet broadening, jet mass, C parameter and two kinds of differ- 
ential two-jet rate. The data cover a large range of the four-momentum transfer Q, which 
is considered to be the relevant energy scale, between 7 GeV and 100 GeV. The Q de- 
pendences of the mean values are compared with second order calculations of perturbative 
QCD applying power law corrections proportional to 1/Q P to account for hadronization 
effects. The concept of power corrections is investigated by fitting simultaneously a non- 
perturbative parameter a p _i and the strong coupling constant a s . 
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1 Introduction 



Hadronic final states in deep-inelastic ep scattering (DIS) ep — > eX offer an interesting en- 
vironment to study non-perturbative hadronization phenomena and predictions of perturbative 
QCD over a wide range of momentum transfer Q in a single experiment. A major limitation 
comes from the treatment of the hadronization of partons, usually modelled by phenomeno- 
logical event generators. Recent theoretical developments suggest that these non-perturbative 
hadronization contributions may be described by 0(1/Q P ) power law corrections P, |2p with 
perturbatively calculable coefficients relating their relative magnitudes. Fragmentation models 
are not needed. 

First results on an analysis of mean event shape variables as a function of Q in terms of 
power corrections and the strong coupling constant a s (Mz) have been published by the HI 
collaboration 0]. It could be shown that power corrections can be successfully applied to the 
variables thrust and jet mass, but they failed to describe the observed jet broadening. In order 
to further test the concept of power corrections, the previous work is considerably extended in 
the present paper. The analyzed integrated luminosity at high Q is tripled, the data correction 
methods are refined and additional event shape variables are investigated. Theoretical progress 
has come from calculations of two-loop effects and the problem of jet broadening has been 
revisited. A comprehensive study of power corrections to the mean values of the event shape 
variables thrust, jet broadening, jet mass, C parameter and differential two-jet rates will be 
presented. The data cover a large kinematical phase space of 7GeV < Q < 100 GeV in 
momentum transfer and 0.05 < y < 0.8 in the inelasticity y. 



2 Event Shapes 

2.1 The Breit Frame 

Event shape analyses in deep-inelastic scattering are based on the observation of the hard scat- 
tering of a quark or gluon which has to be isolated from the target (proton remnant) fragmenta- 
tion region. A particularly suitable frame of reference to study the current region with minimal 
contamination from target fragmentation effects is the Breit frame. Consider, for illustration, 
ep scattering in the quark parton model. In the Breit system the purely space-like gauge boson 
j/Z with four-momentum q = {0, 0, 0, — Q} collides with the incoming quark with longitudi- 
nal momentum p™ z — Q/2. The outgoing quark is back-scattered with longitudinal momentum 
P°qz = ~Q/% while the proton fragments in the opposite hemisphere carrying longitudinal mo- 
mentum p rz — Q/2 ■ (1 — x)/x, where x is the fractional momentum of the struck quark in the 
proton with momentum P. Employing the boson direction as a 'natural' axis, the boost into the 
Breit frame, defined by 2xP + q = 0, thus provides a clean separation into a current and a target 
hemisphere and may be used to classify event topologies at a scale Q. Higher order processes 
like QCD Compton scattering and boson gluon fusion modify this simple picture. However, 
without knowledge of the detailed structure of the hadronic final state, e.g. jet directions, the 
proton remnant is maximally separated from the current fragmentation. 
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The kinematic quantities needed to perform the Breit frame transformation are calculated 
from the scattered lepton (E e >,O e i) and hadron measurements (E^Oi) where i runs over all 
hadronic objects:[| 

Q 2 = 4E e E e , cos 2 ^ , (1) 

y = y e = l-^sin 2 ^ for y e >0.15, (2) 

E e 1 

E^i(l-cos^) 

V = Vh = = for y e <0.15, (3) 

with E e = 27.5 GeV and E p = 820 GeV being the beam energies. The inelasticity y = y e 
is chosen for sufficiently large values. However, since the resolution in y e degrades severely at 
low values, y = y h is taken if y e < 0.15. This procedure ensures least uncertainty in the Lorentz 
transformation to the Breit frame. 



2.2 Definition of Event Shape Variables 

The dimensionless event shape variables thrust, jet broadening, C parameter and jet mass are 
studied in the current hemisphere (CH) of the Breit system. The sums extend over all particles h 
of the hadronic final state in the CH with fo ur- momenta ph = {Eh, Ph}- The current hemisphere 
axis n = {0, 0, —1} coincides with the boson direction. The following collinear- and infrared- 
safe definitions of event shape variables^ are used: 

Thrust t = 1 — T measures the longitudinal momentum components projected onto the 
current hemisphere axis 



E |pfe- n l E \p* 



h 



r = 1 ^ j— = 1 r • (4) 

L \Ph\ L ml 



tGCH heCH 



Thrust re = 1 — Tq uses the direction nj- which maximizes the sum of the longitudinal 
momenta of all particles in the current hemisphere along this axis 

E |Ph ■ n 'l E \Ph ■ n T | 
t c = 1 - max — =-: — — = 1 =-: — : — . (5) 

n'.n' 2 = l E I Phi E IPfcl 

heCH hecn 

This definition is analogous to that used in e + e _ experiments and represents a mixture of lon- 
gitudinal and transverse momenta with respect to the boson axis. 



'Polar angles 8 are defined with respect to the incident proton direction. 

2 Note: The notation of event shape variables is different from the previous analysis [^J, but more transparent. 
All indices are dropped except for tq- The normalization is always performed with respect to the sum of momenta 
or the total energy in the current hemisphere and not Q as was done previously for p. 
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The Jet Broadening B measures the scalar sum of transverse momenta with respect to the 
current hemisphere axis 



E I Pa x n l E IP-l 



h\ 



B = o i — F = o i — r (°) 

2 E |Phl 2 E |p/il 



heCH 



The squared Jet Mass p is normalized to four times the squared total energy in the current 
hemisphere 



( E p h f 

hgCH 

( 2 E E h y • 



(7) 



The C Parameter is derived from the eigenvalues Aj of the linearized momentum tensor jfe 

e piph/\ph\ 



hecn 



E Ip 



h 



heCH 

C7 = 3(A 1 A 2 + A 2 A 3 + A 3 A 1 ). (8) 

The real symmetric matrix jfc has eigenvalues Aj with < A 3 < A 2 < Ai < 1. It describes an 
ellipsoid with orthogonal axes named minor, semi-major and major corresponding to the three 
eigenvalues. The major axis is similar but not identical to nr. If all momenta are collinear then 
two eigenvalues and hence C are equal to zero. 

Higher order processes may lead to event configurations where the partons are scattered 
into the target hemisphere and the current hemisphere may be completely empty except for 
migrations due to hadronization fragments. In order to be insensitive to such effects and to keep 
the event shape variables infrared safe [Q] the total available energy in the current hemisphere 
has to exceed 20% of the value expected in the quark parton model 

Ecu = J2 Eh > ^/ 10 • (9) 

heCH 

Otherwise the event is ignored. This cut-off is part of the event shape definitions, its precise 
value is not critical. 

The event shapes defined in the current hemisphere may be distinguished according to the 
event axis used. Thrust r and the jet broadening B employ momentum vectors projected onto 
the boson direction. Thrust tq and the C parameter calculate their own axis, while the jet mass 
p does not depend on any event orientation. 

Another class of event shapes investigates the number of (n+1) jets found in an event, where 
+ 1 denotes the proton remnant. Jets are searched for in the complete accessible phase space, 
i.e. in both the current and target hemispheres of the Breit frame. Two schemes of jet definitions 
are applied: the Durham or k t algorithm and a factorizable Jade algorithm [||] adapted to 
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DIS. Both jet finding procedures introduce two distance measures: one for distances between 
two four- vectors, y^, and another one for the separation of each particle from the remnant, y, 
The following distance measures are used: 



IT • 



Durham or k t algorithm 



2min(E 2 ,Ej)(l-cos%) 
VU = , (10) 



factorizable Jade algorithm 



2^ 2 (l-cos^ r ) 

y* = -Q2 > ( n ) 



2^(1-CQS%) 
Vij — q 2 , UA> 



2EixE p (l - cos9 ir ) 

Vir = q- 2 , (13) 

where 9 m is the angle between the two momentum vectors. Since the direction of the proton 
remnant coincides with the +z axis for the coordinate system chosen here 9 ir simplifies to the 
polar angle 9{. The pair with the minimal yy or y ir value of all possible combinations is selected 
to either form a new pseudo-particle vector or to assign the particle % to the remnant. The whole 
procedure is repeated until a certain number of jets is found. The event shape variables y^ t (k t 
algorithm) and yfj (factorizable Jade algorithm) are defined as that y value or y ir where the 
transition from (2 + 1) — > (1 + 1) jets occurs. 

Throughout the paper the symbol F will be used as a generic name for any event shape 
variable defined above. Note that for all of them F — > in case of quark parton model like 
reactions. Theoretical calculations of event shape distributions and means will be discussed in 
section |5[ 



3 HI Detector and Event Selection 



3.1 The HI Detector 

Deep-inelastic ep scattering events were collected during the years 1994 — 1997 with the HI 
detector [Q] at HERA. Electrons or positrons with E e = 27.5 GeV collide with E p = 820 GeV 
protons at a center of mass energy of y/s = 300 GeV. Only calorimetric information is used 
to reconstruct the final state. The direction of the scattered lepton and the event vertex are 
obtained by exploiting additional information from the tracking detectors. The calorimeters 
cover the polar angles 4° < 9 < 176° and the full azimuth. 

The calorimeter system consists of a liquid argon (LAr) calorimeter, a backward calorimeter 
and a tail catcher (instrumented iron yoke). The LAr sampling calorimeter (4° < 9 < 154°) 
consists of a lead/argon electromagnetic section and a stainless steel/argon section for the mea- 
surement of hadronic energy. A detailed in situ calibration provides the accurate energy scales. 
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The lepton energy uncertainty in the LAr calorimeter varies between 1% in the backward region 
and 3% in the forward region. The systematic uncertainty of the hadronic energy amounts to 
4%. A lead/scintillator electromagnetic backward calorimeter (BEMC) extends the coverage at 
large angles (155° < 9 < 176°) and is used to measure the lepton at Q < 10 GeV with a preci- 
sion of 1% for the absolute calibration. Since 1995 the backward region has been equipped with 
a lead/scintillating fibre calorimeter improving the uncertainty in the measurement of hadronic 
energies in the backward region from 15% to 7%. The instrumented iron flux return yoke is 
used to measure the leakage of hadronic showers. 

Located inside the calorimeters is a tracking system which consists of central drift and pro- 
portional chambers (25° < 9 < 155°), a forward track detector (7° < 9 < 25°) and a backward 
proportional chamber (155° < 9 < 175°). In 1995 the latter was replaced by backward drift 
chambers. The direction of the scattered lepton is determined by associating tracking informa- 
tion with the corresponding electromagnetic cluster. The lepton scattering angle is known to 
within 3 mrad. The tracking chambers and calorimeters are surrounded by a superconducting 
solenoid providing a uniform field of 1.15 T inside the tracking volume. 



3.2 Event Selection 

The DIS data are divided into a low Q event sample (Q = 7 — 10 GeV, lepton detected in 
BEMC) and a high Q event sample (Q = 14 — 100 GeV, lepton detected in LAr calorimeter) 
which in turn are subdivided further into eight bins in Q: 7 — 8 GeV, 8 — 10 GeV, 14 — 16 GeV, 
16-20 GeV, 20 - 30 GeV, 30 - 50 GeV, 50 - 70 GeV and 70-100 GeV. The following event 
selection criteria ensure a good measurement of the final state and a clean data sample: 

1. The energy of the isolated scattered lepton has to exceed E e i > 14 GeV within 157° < 
9 e , < 173° for the low Q sample and E e , > 11 GeV within 30° < e , < 150° for the high 
Q sample respectively. The calorimetric lepton trigger efficiencies are above 99% [[7| §]. 

2. The inelasticity y is well measured by requiring 0.05 < y e < 0.8 (using the lepton) and 
0.05 < yh (using the hadronic energy flow). This criterion suppresses photoproduction 
events with a misidentified lepton. 

3. The 'quark' direction as calculated from the scattered lepton in the quark parton model 
corresponds to the — z axis of the Breit frame. A minimal value of 9 q > 20° in the labo- 
ratory system ensures a sufficient detector resolution in polar angle after transformation 
into the Breit frame. 



4. A minimal energy in the Breit current hemisphere of Q/10 (see eq. (P|) of section 12^) is 



essential to keep the event shapes r, B, tq, p and C infrared safe. This is part of their 
definition. 

5. To avoid unphysical peaks at zero for rc, p and C at least two hadronic objects are re- 
quired. Events containing only one such object are not quark parton model like but are 
due to leakage into the current hemisphere. 
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6. The total energy in the forward region (4° < 6 < 15°) has to be larger than 0.5 GeV 
to reduce the proportion of diffractive events which are not included into the theoretical 
description of the data. 

7. Hadron clusters have to be contained in the calorimeter acceptance of 5.7° < 6 h < 170° 
avoiding the edges close to the beam pipe. The hadronic energy measured in the backward 
region 6h > 170° has to be less than 10 GeV in order to exclude poor measurements. 

8. The total longitudinal energy balance must satisfy 30 GeV < Ei(l— cos 0j) < 65 GeV 
in order to suppress initial state photon radiation and to further reduce photoproduction 
background. 

9. The total transverse momentum has to be \p±\ < 7.5 GeV (low Q sample) and \p±\ < 
15 GeV (high Q sample), respectively, in order to exclude badly measured events. 

10. The energy measurement of the lepton has to be consistent with that derived from the 
double angle method [Q] \(E e > — E da )/ E da \ < 0.25 in order to further suppress events 
strongly affected by QED radiation. E da is calculated from the directions of the lepton 
and the hadronic energy flow. 

11. An event vertex has to exist within 3a of the nominal z position of the interaction point 

\z v — (z v ) | < 35 cm. 

12. Leptons pointing to dead regions of the LAr calorimeter, i.e. ±2° around 0-cracks be- 
tween modules or ±5 cm around ^-cracks between wheels, are rejected in order to ensure 
a reliable measurement. 

The event selection criteria can be separated into phase space cuts, nos. 1—4, representing 
the common requirements for data and theory, and data quality cuts, nos. 5 — 12. Note that the 
cuts nos. 5 and 6 are always applied except for the perturbative QCD calculations described in 
section |5] where they do not make sense. Depending on the theoretical model to compare with 
they may be considered as phase space cuts as well. Not all cuts affect both the low and high 
Q data samples, as can be seen from the distribution of events in the x — Q 2 plane shown in 
figure [I]. 

The contamination from photoproduction background is estimated to be less than 3% in the 
low Q sample and negligible at higher values of Q. Residual radiative effects are accounted for 
by the data correction procedure described in section |J 

The final data samples consist of 9, 761 events at Q = 7 — 10 GeV taken in 1994 with an 
integrated luminosity of £ = 3.2 pb -1 and 42, 607 events at Q = 14 — 100 GeV corresponding 
to C = 38.2 pb -1 recorded from 1994 - 1997. 



4 Measurement of Event Shapes 

The aim of the analysis is to present measured event shape variables and to compare with sec- 
ond order calculations of perturbative QCD (pQCD) supplemented with analytical power-law 
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X 

Figure 1 : Distribution of selected events in the x — Q 2 plane. Only a small part of the data 
corresponding to a luminosity of 3.2 ptT 1 is shown. The curves indicate the phase space cuts 
in E e >,9 e >,y and 9 q . 

corrections to account for hadronization effects. The experimental observables need to be cor- 
rected for detector effects and QED radiation. Such a procedure involves Monte Carlo event 
generators and hence introduces some model dependencies which are taken into account in the 
analysis of systematic uncertainties. They are much smaller than those obtained in approaches 
which try to unfold from the data to a 'partonic final state' which is directly comparable to 
pQCD. 



4.1 Event Simulation 

A detailed detector simulation is done by using event samples generated with the DjANGO JT0| ] 
Monte Carlo program. The modelling of parton showers involves a colour dipole model as im- 
plemented in Ariadne []TT|]. Alternatively, Django in combination with Lepto [|12|] without 



soft colour interactions has been used to investigate the model dependence!] The hadronization 
of partons uses the Jetset [JT3J] string fragmentation. QED radiation on the lepton side, includ- 
ing real photons and virtual 1-loop corrections, is treated by Heracles [[14]]. Both DIS pack- 
ages provide a reasonable description of the measured event shape distributions and are used to 



correct the data. The event generators Lepto with soft colour interactions and Herwig [ ]T5| ] 



including soft colour interactions in the Lepto version employed spoils the description of the HI data. 
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as stand-alone programs without QED radiation differ considerably in their predictions for the 
hadronic final state and serve to cross-check the unfolding methods. 

4.2 Data Correction Procedure 

The data analysis proceeds in two steps. First the data are corrected for detector effects within 
the phase space described in section [^2] applying different techniques. In a second step QED 
radiation and acceptance corrections due to the beam hole (cut no. 7) are taken into account. 

The most reliable method considered to correct event shape distributions for detector effects 
is found to be a Bayesian unfolding procedure [p)|]. This technique exploits Bayes' theorem on 
conditional probabilities to extract information on the underlying distribution from the observed 
distribution. Although some a priori knowledge on the initial distribution is required — it may 
even be assumed to be uniform in the case of complete ignorance — the iterative procedure is 
very robust and converges to stable results within three steps. The program takes correlations 
properly into account. 

Alternatively, the matrix method employed in the previous publication |J3p and simple cor- 
rection factors applied either bin-by-bin to the distributions or directly to the mean values have 
been used. They serve to estimate systematic effects. The performance of these correction tech- 
niques is checked by studying the spectra and mean values of the event shapes when unfolding 
one Monte Carlo with another and vice versa for various combinations of event generators. In 
general, the Bayes method gives the best results; for details see [[TT]]. 

The remaining corrections account for QED radiation effects and beam hole losses. They are 
applied on a bin-by-bin basis. Non-radiative events are generated using DjANGO with exactly 
the same conditions as before except for the radiative effects being switched off. A detector 
simulation is not required. The corrections are based on the predictions for the hadronic final 
state within the kinematic phase space. 

The derived bin-to-bin correction factors of the event shape spectra are close to one except 
for r, B and y kt . For r and B, which are defined with respect to the boson axis, radiative effects 
are important and non-negligible. The differential two-jet rate y kt is the only variable which is 
sensitive to acceptance losses, particularly at low Q. All other event shapes are almost unaf- 
fected by the beam hole cut as expected for the variables defined in the current hemisphere. Both 
Monte Carlo samples, i.e. Django/Ariadne and Django/Lepto, give consistent results. 

4.3 Results on Event Shape Measurements 

The data correction is performed with the Ariadne event generator as implemented in the 
DjANGO Monte Carlo program. The final results are based on the differential distributions ob- 
tained with the Bayes unfolding and subsequent bin-to-bin radiative correction from which the 
mean values are calculated. The total experimental uncertainties of both the spectra and mean 
values are evaluated in the same way. The statistical uncertainties include data as well as Monte 
Carlo statistics. The systematic uncertainties are estimated by comparing the bin contents and 
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mean of the 'standard' distribution with the values obtained under different conditions or as- 
sumptions. The effects of various correction procedures and the knowledge of the calorimeter 
energy scales are considered. Systematics from the model dependence of the Monte Carlo sim- 
ulation are negligible for the event shapes defined in the current hemisphere alone. In case of 
the y variables the results achieved with Django/Ariadne and Django/Lepto are averaged 
and half of the spread is taken as an estimate of the uncertainty. 

Unfolding uncertainties, being in general asymmetric, are estimated to be half the maximal 
deviation to larger and smaller values respectively due to the alternative unfolding procedures. 
The sum of upper and lower deviation corresponds to half the total spread and so is somewhat 
smaller than twice the standard deviation of a uniform distribution. The influence of the en- 
ergy scale uncertainties is taken into account by repeating the whole analysis and scaling both 
the electromagnetic and hadronic energies separately upwards and downwards by the appropri- 
ate amount. The discrepancies with respect to the central values are attributed to two further 
asymmetric systematic uncertainties. All three (four in case of y) error sources, i.e. the unfold- 
ing bias, the two energy scales and the model dependence, added in quadrature yield the total 
systematics. 

For the event shapes r, B, r c , p and C, the lepton energy uncertainty, which directly affects 
the boost into the Breit frame, is the largest individual contribution, followed by unfolding 
effects. This can be understood because hadronic systematics cancel between numerator and 
denominator for these variables. The situation is reversed for yjj and y^ t . Here, no cancellation 
occurs and the systematic uncertainty due to hadronic energies is the larger of the two energy 
scale uncertainties. 

The corrected event shape distributions are shown in figures f| and |3] over a wide range of 
(Q) — 7.5 — 81.3 GeV. Although the shapes of the spectra for each variable F are quite dif- 
ferent, their common feature is that they all develop to narrower distributions shifted towards 
lower values of F as the available energy Q increases. It demonstrates that the events become 
more collimated. A comparison with the predictions of pQCD reveals serious discrepancies es- 
pecially at low Q leading to the necessity to include hadronization corrections. Such corrections 
are discussed in the next section for the mean values of the distributions. 

The measured mean values are listed in table [[]. Note that the quoted uncertainties do not 
give the correlations due to common systematic error sources. Such correlations are, however, 
properly taken into account in the QCD analysis (see section 



5 Theoretical Framework 

Infrared and collinear safe event shape variables are presently calculable in DIS up to next- 
to-leading order QCD. Several programs are available. However, for comparisons to real ex- 
perimental situations this is insufficient and a phenomenological description of hadronization 
is needed. Within the concept of power corrections one assumes that a parameterization of 
the leading corrections to the perturbative prediction can be obtained without modelling all the 
details of the hadronization. This leads to the notion of 'universal' power corrections with a 
definite Q dependence, typically 0(1/Q P ), given in analytic form with a calculable coefficient 
for each event shape observable. The hope is that such a simplifying approach gives useful 
insight in the interplay of perturbative and non-perturbative effects. 
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Figure 2: Normalized differential distributions of the event shapes r, B and 
t c . HI data (symbols) are compared with Disent NLO calculations (curves) us- 
ing the MRSA' parton density functions with a s (M z ) = 0.115. The error bars 
represent statistical and systematic uncertainties. The spectra given at (Q) = 
7.5 GeV, 8.7 GeV, 15.0 GeV, 17.8 GeV, 23.6 GeV, 36.7 GeV, 57.7 GeV and 81.3 GeV 
(from top to bottom) are multiplied by factors of 10 n (n = 7, . . . ,0). 
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Figure 3: Normalized differential distributions of the event shapes p, C, yfj and 
Dk t - HI data (symbols, see fig. are compared with Disent NLO calculations 
(curves) using the MRSA' parton density functions with a s (Mz) = 0.115. The er- 
ror bars represent statistical and systematic uncertainties. The spectra given at (Q) = 
7.5 GeV, 8.7 GeV, 15.0 GeV, 17.8 GeV, 23.6 GeV, 36.7 GeV, 57.7 GeV and 81.3 GeV 
(from top to bottom) are multiplied by factors of 10 n (n = 7, . . . , 0). 
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(Q) 1 GeV 


(r) 


(B) 


7.5 


n A A no _i_ n nnoo +0.0111 

0.4402 ± 0.0082 I0.0122 


n opo/i _i_ n nno/i +0.0046 

0.3624 ± 0.0034 _o.oo38 


o. f 


n ac\~\ i _i_ n nnnn +0.0196 
U.4U1/ ± U.UUyU _o.oo80 


n Q/Iqc; _i_ n nn/in +0.0076 
U.o4o5 ± U.UU4U _o.oo37 


15.0 


0.3052 ± 0.0034 ±8:88?5 


0.2921 ± 0.0017 ±8:8833 


17.8 


0.2762 ± 0.0029 ±8:8859 


0.2760 d= 0.0016 ±8:8o34 


23.6 


0.2279 ± 0.0031 ±8;So?i 


0.2452 ± 0.0018 ±8:88H 


36.7 


0.1814 ±0.0049 ±8:SS8s 


0.2094 ± 0.0031 ±8:8853 


57.9 


0.1330 ± 0.0089 ±8:8892 


0.1717 ±0.0062 ±8:8ooe 


81.3 


0.0984 ±0.0130 ±8;88tf 


0.1346 ±0.0088 t° mll 



(Q) 1 GeV 


(rc) 


(P) 




7.5 


0.1637 ±0.0032 ±8looio 


o.iii5 =t 0.0019 ±8:88^3 


0.5601 ± 0.0082 ±8:8873 


8.7 


0.1600 =b 0.0037 ±o:8oii 


0.1044 ± 0.0021 ±8:88oo 


0.5524 ± 0.0094 ±8:8854 


15.0 


0.1333 ± 0.0013 ±C;goii 


0.0872 ± o.ooo7 ±8:8821 


0.4824 ± 0.0034 ±8:885? 


17.8 


0.1263 ±0.0011 ±8:8824 


0.0826 ± 0.0007 ±°;{™i? 


0.4621 ± 0.0030 ±8:88eg 


23.6 


0.1098 ± 0.0012 loiooie 


0.0714 ± 0.0007 laooie 


0.4112 ± 0.0033 ±8:8875 


36.7 


0.0985 ±0.0021 ±8:8823 


0.0634 ± 0.0012 t° S 


0.3644 ± o.oo58 ±8:881 


57.9 


0.0834 ± 0.0040 ±8loo4e 


0.0518 ± 0.0023 ±8:882! 


0.3127 ±0.0122 ±8:8?3i 


81.3 


0.0663 ± 0.0057 ±£8825 


0.0410 ± 0.0034 ±8:88?e 


0.2529 ± 0.0173 ±8:8088 



(Q) 1 GeV 


(Vfj) 


(Vk t ) 


7.5 


0.1598 ±0.0026 ±8:8849 


0.3088 ± 0.0065 ±8:8?4 5 


8.7 


0.1497 ±0.0030 ±8:8849 


0.2320 ± 0.0062 ±8:889? 


15.0 


0.1260 ± 0.0015 ±8:881 


0.1349 ±0.0025 ±8:88es 


17.8 


0.1180 ±0.0014 ±8:881 


0.1147 ±0.0021 ±8:8859 


23.6 


0.1049 ±0.0015 ±8:8852 


0.0940 ± 0.0022 ±8:8844 


36.7 


0.0861 ± 0.0025 ±8:8849 


0.0627 ± 0.0028 ±8:8822 


57.9 


0.0785 ± 0.0053 ±8:88ee 


0.0463 ± 0.0044 ±8:881 


81.3 


0.0608 ± 0.0071 ±8:8857 
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Table 1: Corrected mean values of the event shapes as a function of Q. The first uncertainty is 
statistical, the second systematic. 
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Figure 4: Mean values of r (left) and C (right) versus Q in four different bins of x calculated 
with Disent. The lines connect the means belonging to the same x bin. 



5.1 pQCD Calculations 



The mean value of an event shape variable F can be written in second order perturbative QCD 

as 



pert 



ciO, Q) a s (fi R ) + 



c 2 (x, Q) + ^ In ^ci(x, Q) 



(14) 



where ji R is the renormalization scale, (5q = 11 — 2/3 Nf and Nf = 5 is the number of active 
flavours. In contrast to e + e~ annihilation, where the coefficients ci and c 2 are constant, in 
DIS they depend on x according to the parton density functions and the accessible x-range 
at different values of Q (see figure [1]). Examples for a strong and a weak variation with x 
of the mean values (t) and (C) versus Q are presented in figure |]. Average values of the 
coefficients c\ and c 2 are calculated separately for every Q bin with its specific range of re- 
values. This approximation by step functions is the origin of the steps exhibited by the curves 
in figures |5] and ^. The slope within one bin is due to the variation of a s with Q. In the 
previous publication [[3J], average coefficients determined for the complete range in Q have been 
employed. 



The event shape means are evaluated at the scale fi R = Q with the Disent [ |18fl program 
which treats deep-inelastic ep scattering to 0(a 2 s ) in the MS scheme and employs the subtrac- 
tion method for the necessary integrations. In order to test the reliability of the perturbative 
predictions, careful comparisons of Mepjet [jI9p, Disent and Disaster++ [[ZUp have been 



carried out [plQ. Due to intrinsic restrictions of the integration technique applied in Mepjet 
(phase space slicing), however, it can not be used for the mean event shapes. In general, good 
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agreement at the percent level is observed but some discrepancies have been revealed. These 
are partially understood and an improved Disent version has been employed. The effects are 
hardly visible in the event shape spectra, but the mean values have increased considerably in 
the low Q region compared to the calculations used in [Q]. The worst case is the jet broadening 
B whose mean value at Q = 7.5 GeV rose by about 14%, being still slightly different from the 
DlSASTER+4- calculation.^ 



The parton density distributions MRSA' [E2fl with a s (M z ) = 0.115 are used as standard. 



Other sets, i.e. MRST99 [ p3Q and CTEQ4 Q24[], are investigated as well to estimate systematic 
uncertainties. 



5.2 Power Corrections 

Hadronization effects on event shapes are treated within the concept of power corrections [[l],[2|]. 
The observable mean values can be written as 

(F) = (F) port + (F) pow , (15) 

with (F) per given by eq. (|T4|). The hadronization contributions (F) pow are expected to be 
proportional to 1/Q P with exponents p = 1 or p = 2 depending on the observable. Two types 
of parameterizations for (F) pow will be investigated. 

In a simplistic approach inspired by the longitudinal phase space or tube model, described 
e.g. in [|25|], one has 

(F) pow = ^f or (F) pow = ^f, (16) 

where X^p and A2,f are constants. One expects Ai to be in the order of ~ 1 GeV and terms 
with exponents p > 1 to be negligible except for the differential two-jet rate y^ t . Even within 
this simple model it is possible to derive approximate relations between the constants Ai for 



different event shapes, e.g. Ai T ~ 2 Ai i|0 , consistent with e + e data [25]. 



In the model pioneered by Dokshitzer and Webber [[I], the idea is to attribute l/Q p power- 
law corrections to soft gluon phenomena associated with the behaviour of the strong coupling 
at small scales. This leads to the notion of a universal infrared- finite effective coupling a e g(/j, R ) 
which replaces the perturbative form in the infrared region fi R < [i 1 where fj, z , the infrared 
matching scale, has to fulfil Aq CD <C /i 7 ~ 2 GeV <C Q. At the expense of one new non- 
perturbative parameter « p _i(/i / ), corresponding to the (p — l)th moment of the effective cou- 
pling a c ff(fi R ) when integrated from up to the power corrections to all event shapes with 
the same exponent p can be related via 



(F) pow = a F V, (17) 



v = ACf m> 




a p -i 


7T p 


K^r) 





lOx) - a s (fx R ) - ^- ( In ^ + ^ + - ) a 2 s 
27T V (J p ' 



(18) 



4 Artificially increasing the DlSENT predictions for the two low Q means of the jet broadening by 10% and 5% 



respectively improves the consistency of the fits described in section 5.2, especially with respect to a s . 
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where Cp = 4/3, K = 67/6 — n 2 /2 — 5/9 Nf and Nf = 5 as in the perturbative part. Again, 
fi R is identified with Q in this study. The subtractions proportional to a s and a 2 s serve to avoid 
double counting. 

The coefficients dp, given in table depend on the observable F but can in principle be 
derived from a perturbative ansatz, although not yet available for the variable . Some of them 
have changed considerably — e.g. a# and ac by factors of 4 and 2 respectively — from their 
original values [ |2~8l ] applied in previous investigations [[3], ^9). Ambiguities in calculating the 
ap predictions could be resolved with the advent of a two-loop analysis provided an additional 



common coefficient, the Milan factor M. [26, 27] with 



2 . , 2 / 4.725 - 0.104 N f \ 
M = -M = -il + - L ) (19) 

is applied. The numerical value is M! ~ 0.95 for Nf = 3 flavours relevant for gluon radiation 
at low scales. [] 

For all event shapes under study a power suppression exponent of p = 1 is expected except 
for the two-jet rate ?/fc t where p = 2. The power correction parameter «o is estimated to be 
~ 0.5 whilst «i is essentially unknown. 

The jet broadening B is expected to behave differently from the other event shapes and 
eq. ( |T8| ) should contain an additional enhancement. The originally proposed factor of ln(Q / Q ) 
[jZSp, with Q ~ CK/U/) an unknown scale, could not be supported by the HI data [|3|]. This 
observation stimulated a theoretical reexamination leading to the following power correction 
for the jet broadening [TJO] 



(B>P "' = a - a ^- B ' 2v /^w 1 I^).^) + i-^ + "° |P - (20) 



3 

where a B = 1/2, ?7 = —0.614 and a s has to be evaluated at the scale • fi R . The en- 
hancement term a' B is substantial with a slow variation of 1.6 — 2.2 over the measured Q range. 
Strictly, this formula has been derived for e + e~ annihilation, but should be applicable for DIS 
as well. The accuracy of the coefficient, however, is only of the order of 1 [[$0|] ! All coeffi- 
cients dp are considered to be input parameters and systematic uncertainties do not account for 
approximations, e.g. the neglect of quark mass effects, inherent in their derivation. 



6 QCD Analysis of Event Shape Means 

In order to get an impression of the impact of hadronization, figures [2] and ^] show the ex- 
perimental event shape spectra in comparison with calculations of perturbative QCD. At high 
values of Q the effects are small, as expected, and the unfolded and partonic spectra approach 
each other. At low Q, data and calculations look very different and non-perturbative effects 

5 Note that the original derivation, which lead to M 1 ~ 1.14, has recently been corrected [^7[. 
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become prominent. It is particularly interesting to note that both two-jet rates yjj and y^ t al- 
ready exhibit small hadronization corrections at modest momentum transfer, a characteristic 
very different from the other variables. 

In this section a comprehensive study of event shape means will be presented in order to 
pin down the analytical form and magnitude of power-law corrections and to test their universal 
nature. The strategy is to investigate the Q dependence of (F) as given in table [I] assuming 
the ansatz of eq. (JT3|) . The perturbative part is given by the QCD expression of eq. (JT4|). Two 
variants of power corrections (F) pow will be tested: the tube model and the approach pioneered 
by Dokshitzer and Webber. 

6.1 l/Q p Fits 

Keeping a s (Mz) fixed to 0.119 and parameterizing hadronization contributions by eq. ( |T6| ) 
yields the fit results for \i t p and \2,f given in table ||. With the exception of r, B and y/j 
acceptable results can not be achieved in the case of the \\/Q term. Note that for r and B Ai 
complies with the expectations. Applying A2/Q 2 corrections the x 2 values worsen dramatically 
and lead to the rejection of this ansatz. 



(F) 


A ljF /GeV 


xVdof 


A 2 , F /GeV 2 


rVdof 


a s (M z ) 


X 2 /dof 


(r) 


0.71 ±0.03 


2 


7.2 ±0.3 


23 


0.131 ±0.001 


1 


(B) 


0.55 ±0.02 


1 


5.2 ±0.2 


28 


0.134 ±0.001 


24 


(re) 


0.73 ±0.01 


29 


7.1 ±0.2 


212 


0.152 ±0.001 


2 


(P) 


0.54 ±0.01 


38 


5.2 ±0.1 


283 


0.160 ±0.001 


1 


(C) 


2.33 ±0.03 


66 


20.8 ±0.4 


327 


0.148 ±0.001 


5 


{Vfj) 


-0.13 ±0.01 


2 


-1.2 ±0.1 


7 


0.113 ±0.001 


1 


(Vkt) 


-0.32 ±0.02 


7 


-3.1 ±0.3 


21 


0.110 ±0.001 


9 



Table 2: Results of fits to event shape means. Left: fit of Ai,i? with a s {Mz) = 0.119, center: fit 
of A2,f with a s (Mz) = 0.119 and right: fit of a s (Mz). Uncertainties are statistical only. 

Allowing for test purposes a variation of a s (Mz) in the perturbative expression while ne- 
glecting power corrections, some of these fits actually work but result in large discrepancies of 
the strong coupling with respect to a world average of a s (M z ) = 0.119 ± 0.002 [plj]. Neither 
a power-like correction according to eq. dl~6l ) nor pure pQCD, exploiting the logarithmic Q de- 
pendence of the strong coupling, are sufficient to describe all data. However, fitting a s (M z ) and 
\i t F or A 2j f simultaneously, does not lead to satisfactory results either. Due to an extreme anti- 
correlation between the two parameters there is a tendency to minimize the power contribution 
at the cost of unphysical shifts of a s (M z ). Only B and y^j produce correlated but reasonable 
numbers. Details can be found in []17[j. These studies suggest that some form of combined 
power-like and logarithmic Q dependence is needed to account for the observed medium (r, 
B), large (t c , p, C) and small {yt t , yjj) hadronization corrections. 
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/ r~i\ 

(F) 


CLf 
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a v -\\\i l = 2 GeV) 
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0.5 


-98 


(B) 


1/2 • a' B 


1 


fl ^7 +0-017 +0.028 
U.OO 1 _ .012 -0.039 


niiio +0.0036 +0.0049 
U.iliO _ .0028 -0.0051 


0.7 


-69 


(P) 


1/2 


1 


D ^07 +0.009 +0.050 
V.OJI _ .010 -0.057 


fl 1 37/1 +0.0024 +0.0110 
U.iCW^ -0.0032 -0.0096 


1.1 


-32 


(re) 


1 


1 


n c;n^ +o.oos +0.043 

U.OUO _ .010 -0.047 


n 1 01 n +0.0023 +0.0098 
U.IOIU -0.0028 -0.0089 


1.2 


-22 


(C) 


3tt/2 


1 


fl AA7 +0-005 +0.032 
V-ltl _o.007 -0.038 


n 1 +0-0016 +0.0103 

U.lOUl -0.0020 -0.0091 


0.8 


+36 


(Vfj) 


1 


1 


0.28 ±8^ 


nine; +0.005 

U.1UO -0.006 


0.8 


-72 


(Vfj) 


-0.2* 


1 


n 07 +0.20 

u '°' -0.21 


n -1 1 fi +0.008 
U.llO -0.009 


0.6 


+98 


(Vk t ) 


It 
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u ' 00 -0.04 


001 +0 022 
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-98 
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It 
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1-50 j£S 


n OQQ +0 007 
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3.6 


-92 


(Vkt) 


-50* 


2 
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U.lZ^t -0.014 


0.6 
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Table 3: Results of fits a la Dokshitzer-Webber for the event shape means. The coefficients 
ap and exponents p of the power corrections are given as well. The first uncertainty contains 
statistics and experimental systematics, the second is an estimate of theoretical uncertainties 
(omitted for y variables), k denotes the correlation coefficient between and a s (M z ). The 
starred coefficients of the y variables are derived from a fit procedure, whereas the coefficients 
marked with t are trials/guesses. 



6.2 Fits in the approach initiated by Dokshitzer and Webber 



Fits in the approach initiated by Dokshitzer and Webber are performed separately for each event 
shape. The data are very well described by pQCD plus these analytical power corrections as 
shown in figures |5| and ELR The fit results for the power correction parameters a p _i and the 
strong coupling a s (M z ) are compiled in table |} For a discussion of the quoted uncertainties 
see section |53| . Theoretical systematics for the fit parameters of the differential two-jet rates 
yfj and y kt are not given due to the essentially unknown coefficients a F in the power correction 
predictions. 

First the fits to event shape variables defined in the current hemisphere will be discussed, 
i.e. not including the two-jet variables yjj and y^. With the exception of r, correlations are 
reduced compared with the tube model and the x 2 values are reasonable. It is interesting to note 



that the new calculations of the power corrections for the jet broadening, eqs. (JTVp and (|20|X 
are now able to describe the data well. Neglecting the enhancement factor gives fit results of 

a = 0.661 iHg and a s (M z ) = 0.1169 ^ t° S with x^/doi = 0.8. 

The fitted parameters are displayed in the a s — a plane of figure [7[ Note that the ex- 
perimental uncertainties, statistics and systematics, are in general smaller than the theoretical 
uncertainties of 5 — 10%. The parameters a scatter around the expectation of a ss 0.5 within 
about 20% and are compatible with the assumption of universality in the Dokshitzer-Webber 
approach. The spread of the strong coupling constant appears to be uncomfortably large and 
one observes a group of higher a s (Mz) values for those variables which do not make use of 



6 The steps are due to the x dependence of the pQCD calculations (see section 5.1) 
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Figure 5: Mean values (full symbols) of r, B, t c , p and C as a function of Q. The error 
bars represent statistical and systematic uncertainties. The full line corresponds to a power 
correction fit according to the Dokshitzer-Webber approach. The dashed line shows the pQCD 
contribution of Disent in these fits. 
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Figure 6: Mean values (full symbols) of yfj and y kt as a function of Q. The error bars represent 
statistical and systematic uncertainties. Upper part: The full line corresponds to a fit of the 
pQCD calculation without power contribution. Lower part: The full line corresponds to a power 
correction fit according to the Dokshitzer-Webber approach. The dashed line shows the pQCD 
contribution of Disent in these fits. 
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Figure 7: Results of the power correction fits to the mean values of r, B, p, tq and C with 
the contours of x 2 («s,«o) — Xmin + 1 an d X 2 ( a s,®o) = Xmin + 4 including statistical and 
experimental systematic uncertainties. 



the boson axis. Possible explanations are missing higher order QCD corrections, expected to be 
different for each event shape variable, and/or incomplete knowledge of the power correction 
coefficients. Allowing the coefficients ap to vary by arbitrary factors of 2 and 1/2 in order to 
study the effect on a Q and a s (Mz), one observes shifts of the ellipses of figure [7] approximately 
along the main diagonal. Larger coefficients ap induce smaller values of a and a s (M z ) and 
vice versa. 



As already seen in section O] both two-jet rates y defined in the whole phase space of the 



Breit system exhibit only small hadronization effects, much smaller than e.g. thrust. For the 
two-jet rate yjj based on the Jade algorithm an exponent p — 1 and a coefficient a yfJ = 1 
for the power corrections have been given in [Q] . Formally, one obtains an acceptable fit but 
rather low and unreasonable numbers for a and a s (Mz). In particular any value of «o ^ 0.3 
makes no sense within the concept of this model. This leads to the conclusion that the power 
correction coefficient of ref. []2|], which has not been reinvestigated since, may not be correct. 
Instead, more consistent results can be obtained by either neglecting hadronization corrections 
altogether or by taking a small, negative coefficient of a Vf J = —0.2 (see figure || and table |3|). 
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In contrast to the other event shapes the hadronization correction to (y kt ) is expected to 
decrease with a larger power of p = 2 []2|] . In fact a power law behaviour with p = 1 and a Vk = 1 
is strongly disfavoured by the data. Other choices of a ykt can not accommodate for reasonable 
values of ao > 0.3. For a power exponent of p = 2 the coefficient a Vkt is basically unknown 
and has to be determined in addition to oi\. One can try to fit simultaneously ai, a s (Mz) and 
a yk . This set of parameters, however, is strongly correlated. Nevertheless, the fit converges 
properly giving «! = 0.34±0.05, a s (M z ) = 0.125 ±0.006 and a ykt = -52±3 with x 2 /dof = 
0.6 and statistical uncertainties only. Reinserting a value of a Vk = —50 gives the entry in 
table |3J. To prove a power correction term oc 1/Q 2 , however, is currently not possible because 
the observed hadronization effects are small. The two-jet rate data y kt are certainly consistent 
with quadratic power law corrections, but more experimental precision and theoretical input is 
needed. Figure ^| shows fits to the mean values without and with hadronization contributions. 

In summary, the experimental observations are very well described within the approach ini- 
tiated by Dokshitzer and Webber. The analytical form of the power correction contributions 
appears to be adequate. Also the magnitude is of the right order for those event shape vari- 
ables where updated calculations exist, supporting the notion of approximate universal power 
corrections. 



6.3 Systematic Uncertainties 

The procedure to estimate the systematic uncertainties in table |3| is to repeat the fits under 
variation of every prominent systematic effect. The discrepancy compared to the standard result 
is attributed to a corresponding uncertainty. In case of deviations in the same direction for the 
variation of one primary source, e.g. an upwards and downwards modification of an energy 
scale, only the larger one is considered for the evaluation of the total uncertainty. The latter is 
derived from all contributions added in quadrature. An exception is the unfolding procedure 



whose influence is estimated as explained in section |4~3| . The following systematic effects are 
investigated: 

• Experimental uncertainties 

1 . Usage of four correction procedures 

2. Variation of the electromagnetic energy scale of the calorimeters by ±(1% — 3%) 

3. Variation of the hadronic energy scale of the LAr calorimeter by ±4% 

• Theoretical uncertainties 

1. Variation of the renormalization scale //| by factors of 2 and 1/2 

2. Variation of the factorization scale fx 2 F by factors of 4 and 1/4 

3. Variation of the infrared matching scale / u / by ±0.5 GeV 

4. Usage of MRST99 parton density functions [ |2"3| ] with larger and smaller gluon con- 
tributions and usage of MRSA' parton density functions p2|] with strong couplings 



a s (Mz) = 0.105 up to 0.130, different from the standard set 
5. Usage of different parton density functions CTEQ4A2 [ j24] ] with similar a s (M z ) 
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The experimental sources are already discussed in section The renormalization scale 



fi R and the factorization scale \i F are arbitrary since in a complete theory the calculations do not 
depend on any specific choice. But in reality one has only an approximate theory yielding resid- 
ual dependences due to neglected higher orders. To avoid the appearance of large logarithms in 
the computations it is recommended to identify the scales with a process relevant scale chosen 
to be Q in the present analysis. To estimate the effect of higher orders it is conventional to vary 
fi 2 R and fi 2 F by an arbitrary factor of 4. In the case of fi 2 R one has to reduce this factor to 2 because 
of the infrared matching condition A QCD <C /x 7 ~ 2 GeV <C p, R . 

The variation of the infrared matching scale by pi r = (2 ± 0.5) GeV follows the original 
proposal [[|]. Note that this affects only the a s (Mz) uncertainty. The parameter a explicitly 
depends on 

The last two points account for uncertainties in the gluon content of the proton and the fact 
that a s (M z ) has implicitly already been used in deriving the parton density functions which 
may bias the computations. The same is true for the choice of a parameterization for the parton 
density functions. Therefore five alternative sets, two each with different gluons or a s (Mz) 
respectively and one with approximately the same a s (Mz) but another parameterization, are 
chosen for a reevaluation of the Disent calculations. 

The contributions of all systematic uncertainties are presented graphically in figures |]and|9| 
for each of the five event shapes r, B, p, tq and C where the coefficients ap are known. Without 
such a prediction for the two-jet event rates yfj and y^ t the study of systematic effects has been 
performed for experimental uncertainties only. 

From figures |8] and |9| it is obvious that r and B have different properties compared to p, 
tq and C. For the latter the systematic uncertainties are dominated by the variation of the 
renormalization scale. In case of r and B there is no prevailing source of systematics. The 
larger influence of experimental uncertainties can probably be related to the explicit reference 
to the boson axis implied in the definitions of r and B. 

Several additional cross-checks concerning systematic effects are performed, (i) Leptons 
pointing to partially inefficient regions in the LAr calorimeter (see section [3^ ) are removed. 
This selectively diminishes the contribution to certain phase space regions and is corrected for. 
It is checked that the actual influence is negligible even without unfolding, (ii) The power 
correction coefficients do not account for phase space constraints imposed on the scattered 
lepton (cut no. 1). Experimentally they are unavoidable, but for testing purposes the pQCD 
calculations are repeated without these cuts. Except for (y kt ) which increases by ~ 8% at low 
Q all other mean values change by less than 2%. (iii) In the power corrections eqs. (|T8|), < JT9| ) 



and ( PD| ) it is not obvious which number of flavours to take. As default Nf = 5 is used for 
the perturbative part and the subtraction terms and Nf = 3 in the coefficients of the power 
corrections. Repeating the fits with Nf = 3 in both cases increases all ao by ~ 0.03. The 
values of a s (M z ) are almost unaffected. 
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Figure 9: Systematic uncertainties of aof^/ = 2 GeV) and a s (M z ) for r c (top), p (middle) and 
the C parameter (bottom). 
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7 Conclusions 



The event shape variables thrust, jet broadening, jet mass, C parameter and two-jet event rates 
are studied in deep-inelastic ep scattering. The mean values exhibit a strong dependence on the 
scale Q, which can be understood in terms of perturbative QCD and non-perturbative hadroniza- 
tion effects decreasing with some power p of Q. The interest is to test whether the leading 
corrections to perturbation theory can be parameterized without any assumptions on the details 
of the hadronization process. Simplistic models of the form const/Q p with p = 1 or 2 fail. 
The concept of power corrections in the approach initiated by Dokshitzer and Webber, which 
predicts the power p as well as the form and magnitude of the hadronization contributions as 
a function of a new parameter provides a much better and satisfactory description of the 
data. 

The event shape variables defined in the current hemisphere of the Breit frame — r, r c , B, 
p, and C — have sizeable power corrections proportional to 1/Q. Two-parametric fits yield for 
the non-perturbative parameter a ~ 0.5 within 20% supporting the notion of universal power 
corrections for the event shapes. The corresponding, correlated values of the strong coupling 
a s (M z ) show a large spread incompatible within the experimental uncertainties. Possible ex- 
planations are missing higher order QCD corrections, expected to be different for each event 
shape variable, and/or incomplete knowledge of the power correction coefficients. 

Both two-jet rates y, defined for both hemispheres of the Breit system, are almost unaffected 
by hadronization effects. For yjj the conjectured large positive contribution with ajj = 1 is 
ruled out by the data, which instead prefer small negative power corrections. For y kt fits with 
p = 1 do not work properly. The data are consistent with the expectation of a power p = 2, 
but unfortunately more quantitative statements can not be made given the current experimental 
precision and the unknown coefficient a,k t . 

The improved and extended analysis of mean event shape variables in deep-inelastic scatter- 
ing supports the concept of power corrections. In order to achieve a better common description 
of the different event shapes and to get a more coherent picture further theoretical progress is 
needed. 

Similar studies of power corrections, often combining data from several experiments to 
cover a sufficiently large range in Q, have been done for the analogous e + e~ event shape vari- 
ables r, B, p and C [[32"!] . Remarkably, these analyses also find universal parameters a ~ 0.5 
within 20%. But depending on the choice of employed data sets correlated values of a and 
a s (M z ) are found which deviate between the different determinations by more than the quoted 
experimental accuracies. The situation appears to be similar to the one in DIS. 

Note: The experimental data of the event shape spectra shown in figures [Z] and [5] may be 
obtained in tabular form from the HEPDATA data base at Durham. 
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